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We analyze a model system of fermions in a harmonic oscillator potential under the influence of 
a fluctuating force generated by a bath of harmonic oscillators. This represents an extension of the 
well-known Caldeira-Leggett model to the case of many fermions. Using the method of bosonization, 
we calculate Green's functions and discuss relaxation and dephasing of a single extra particle added 
above the Fermi sea. We also extend our analysis to a more generic coupling between system and 
bath, that results in complete thermalization of the system. 
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The interaction between a system and its environment 
is an important fundamental issue in quantum mechan- 
ics. It is at the basis of relaxation phenomena (like 
spontaneous emission), is essential for the measurement 
process, and leads to the destruction of interference ef- 
fects ( "decoherence" or "dephasing"). In the theory of 
quantum-dissipative systemsQ, there are only few ex- 
actly solvable models, most notably the Caldeira-Leggett 
modeljl of a single particle coupled to a bath of harmonic 
oscillators. This is the simplest possible model in which 
friction and fluctuations appear. If the particle is free, 
then this model can be used to study the quantum ana- 
logue of Brownian motion. The model remains exactly 
solvable if the particle moves in a parabolic potential (the 
damped quantum harmonic oscillator). 

However, in many solid state applications, we actu- 
ally consider dephasing of an electron inside a Fermi 
sea. It is difficult to apply the insights gained from 
single-particle calculations in such cases, since the Pauli 
principle may play an important role in relaxation pro- 
cesses. There have been comparatively few detailed 
studies of quantum-dissipative many-particle systems. 
Among them we mention a general discussion of dephas- 
ing in a Luttinger liquid Q , a study of fermions coupled 
to independent baths Q , and a formally exact extension 
of the Feynman- Vernon influence functional to fermions 
In other cases, the Pauli principle has been intro- 
duced "by hand" , by keeping only the thermal part of 
the bath spectrum |(|. 

In this Letter, we study a natural extension of the 
Caldeira-Leggett model to a many-fermion case. The 
model consists of a sea of fermions populating the lower 
energy levels of a harmonic oscillator. We are interested 
in the effects that arise when a bath is coupled to this 
system via a fluctuating spatially homogeneous force. In 
contrast to an analogous system of free fermions Q , the 
bath leads to transitions between levels, with strong ef- 
fects of the Pauli principle. This model might also prove 




Figure 1: Left: Fermions in an oscillator, coupled to a dissi- 
pative bath via a fluctuating force F. Right: Approximately 
equivalent model of chiral fermions on a ring, subject to a 
transverse force. 

relevant to the discussion of cold fermionic atoms in a Id 
harmonic trap under the influence of fluctuations 

of the trapping potential. 

We rewrite and solve the Hamiltonian using the 
method of bosonization, for the case of large particle 
numbers. This enables us to evaluate Green's functions 
and to describe relaxation and dephasing of an extra par- 
ticle added above the Fermi sea. Finally, we will extend 
our model to a more generic type of coupling. 

The model - We consider a system of N identical 
fermions (non-interacting and spinlcss) confined in a one- 
dimensional harmonic oscillator potential (see Fig. 
A fluctuating force F leads to a coupling of the form 
F^2,j Xj, yielding, in second quantization: 

oc 

H = w ° n< ^c„ + H B + 

m ^TT(ci +1 c„ + h.c.) (i) 

The operators c n annihilate fermions in the oscillator lev- 
els n. The bath Hamiltonian Hb describes an infinite 
number of harmonic oscillators, and the force F is a sum 
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over the bath normal coordinates Qj . It is characterized 
fully by its power spectrum (^FF^ . The special case 
of an Ohmic bath, used for Quantum Brownian motion 
2], has (PP^ = (riu>/ir)6{uj c - lu)6{lu) at T = 0, where 

rj = mj is the friction coefficient, 7 the damping rate, 
and u! c the cutoff. As the form of the coupling (JIJ is 
not translationally invariant, the frequency ojq contains 
a stabilizing counterterm|lL la flO| . 

Effectively, the force acts only on the center-of-mass 
(cm.) motion of the particles which, for the harmonic 
oscillator, is independent of the relative motion. Thus, 
in principle our model reduces to a single damped har- 
monic oscillator, analyzed in Ref. 0- However, we are 
interested in single- fermion properties, and not in the col- 
lective cm. motion itself. Although the problem can be 
solved exactly via normal modes of the complete set of 
oscillators and antisymmctrizing with respect to fermion 
coordinates, this procedure gets extremely cumbersome. 
Instead, we employ an approximation for large fermion 
numbers N, which also allows an extension to a more 
generic coupling between system and bath. 

Bosonization - For sufficiently large N the lowest lev- 
els are always occupied (at the given interaction strength 
and temperatures), i.e. excitations are confined to the 
region near the Fermi level. Then we may employ the 
method of bosonization, rewriting the energy of fcrmions 
as a sum over boson modes 0. This is possible since 
the energies of the oscillator levels increase linearly with 
quantum number, just as the kinetic energy in the Lut- 
tinger model of interacting electrons in one dimension 
(for recent reviews see H3)- 

We introduce (approximate) boson operators b q = 
l^nLo cliCn+q (q > 1), which destroy particle-hole ex- 
citations. Then, the Hamiltonian given above becomes 
approximately 




ip{x) dx 



(3) 

The coordinate x does not refer to the motion in the 
oscillator. Rather, wc have effectively mapped our prob- 
lem to a chiral Luttinger liquid on a ring with a coupling 
tx ^cos(a;) (x G [0,2tt[), see Fig. [T] (right). Thus, the 
following results also describe relaxation of momentum 
states in that model. Although a generic discussion of 
dissipative Luttinger liquids has been provided in Q , the 
particular questions we are going to study have not been 
analyzed before. 

The operators tjj(x) may be expressed as[ll|: 

i){x) = K\{x)e^ {x) e l ^ x) = KXe^r , (4) 

with 

The "Klein factor" K annihilates a particle, with 
[K, bp] = and K(t) = Kexp(-iu (N - l)i). We have 
X(x) — exp(i(N — l)x)/y/2w and r = exp(—[<p^,0\/2). 
(The exponent in r diverges, so a formal cutoff at high 
q should be introduced, which will drop out in the end 
result) 

Using Eq. 0, we find for the hole-propagator: 

= ^ J 2 \^'*'-™) (ft( x > it )$(x,0)) dxdx'. 

(6) 

The ^-Green's function is given directly in terms of the 
^-correlator, using Eq. Q: 



H « Lo V qb\b q + J F(h +b\)+H B + E^, (2) 

9=1 

which will form the basis of our analysis. Here Efj — 
u>qN(N — l)/2 is the total energy of the iV-fermion non- 
interacting ground state. Eq. J2J reveals that F only 
couples to the lowest boson mode (q= 1), corresponding 
to the cm. motion. The damped motion of the cm. os- 
cillator can be solved exactly, along the lines of Ref. |2J or 

[l2|. providing us with correlators such as (bi(t)b\(0)^ . 

Derivation of Green's functions - In order to find the 
Green's functions, we have to go back from the boson op- 
erators b q to the fermion operators c„, by employing well- 
known finite-size bosonization identities. In our case, we 
first have to introduce auxiliary fermion operators ip{x): 



^ t (x',t)^(a;,0)\ = !L e in F ((x-x')+u t) / e -4{x',t) e 4(x,0) 

(7) 

(with up = N — 1). The expectation value on the 
right-hand side may be evaluated exactly since the 
system-bath coupling is bilinear. This yields exp(E) 
with: 

E = ~ ((0(x',t) 2 ) + (#c,0) 2 )) + (<Kx',i)#E,0)) 

The correlator of is a polynomial in X = exp(ix) and 
X' = exp(ix') (see Eq. (0). Now the double Fourier in- 
tegral in Eq. JSJ) may be evaluated by expanding exp(E) 

as a series in X and X'. We find that ^cjj/(£)c ra \ is the 
coefficient of X n /X' n in the expansion of 
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JE(X.X',t) 



E 



e iuoH (X/X') k , 



(9) 



where the noninteracting exponent has been subtracted 
in SE = E — E(Q\, which thus contains only the corre- 
lator of the damped cm. mode q = 1. Detailed plots 
of the Green's function will be published elsewhere [lfT |. 
Here we provide the result in the weak-coupling approx- 
imation, where we neglect the bath-induced smearing of 
the equilibrium Fermi level and use an exponential de- 
cay for the cm. motion. Both assumptions are summa- 
rized in (bi(t)bi^ = (b{(t)bi\ = and ^(^(O)^ = 
e -iw t-rt/2 j g renormalized cm. frequency, and 
F is the decay rate, with T = N'y for the Ohmic bath). 
Then, the exact Eq. © yields 



(4(<)c n ) 



m=0 



V (ty 



(10) 



where v(t) = exp(— i(u — ujo)t — Tt/2) — 1. 

We find that the hole (particle) propagator does not 
decay to zero in the limit t — > oo, for any n < uf — 1 (n > 
rip + 2), since v (t) — ► — 1. This is in contrast to the naive 
single-particle picture of complete decay for any level n ^ 
Uf, tif + 1 not directly at the Fermi level (i.e. the result 
suggested by the leading order self-energy). Physically, 
adding a hole (particle) creates an excited many-particle 
state which also contains contributions where the cm. 
mode is not excited, and these will not decay, because 
only the cm. mode is damped. A more generic coupling, 
leading to ergodicity, will be discussed at the end of this 
Letter. 

Time- evolution of density matrix - We now turn to 
the two-particle Green's function in order to learn about 
relaxation of level populations and dephasing. Con- 
sider placing an electron in a superposition of levels 
above the Fermi sea, creating the many-particle state 
Sn ^no^no \FS) at time 0. We assume the levels n 
to be unoccupied. This will hold for n > rip in the 
weak-coupling limit, for which the following results have 
been evaluated. The reduced single-particle density ma- 
trix evolves according to: 



Pnn' (t) 



n ,n' 



* no * r. 



•»' D cI'(*)£n(*)4 ( 



(11) 



We may rewrite the Green's function in Eq. in terms 
of tp(x) (Eq. J3J), leading to a four- fold Fourier integral, 
analogous to Eq. ©. Using Eq. Q), this may be evalu- 
ated by a series expansion in four exponentials exp(ix^ '), 
exp(i?/ )), similar to Eq. ©. We omit the lengthy gen- 
eral formula 0, but dfscuTs a limiting case below. 
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Figure 2: Time-evolution of the single-particle density ma- 
trix p nn '(t) (Eq. mil l, after placing an extra particle in a 
superposition of two states. The radius of each circle gives 
|p nn /(t)|. Levels n,n' = % — 1, . . . ,np + 6 are indicated by 
grid lines. For this example u>o — cj'o = 2T. 



In Fig. |2 we have plotted the resulting time-evolution 
of the density matrix for the case of an equal superposi- 
tion of two levels, = \I>„ 2 = 1/V2, at T = 0. The 
population of the highest occupied states in the Fermi sea 
(lower left of panels) decreases, because these fermions 
become partly excited at the expense of the extra parti- 
cle, due to the effective interaction mediated by the bath. 
Moreover, the particle does not decay all the way down 
to the lowest unoccupied state np + 1. Rather, in the 
long-time limit, the excitation is distributed over a range 
of levels above the Fermi level, up to the initial levels 
711,2- Again, this is because only the cm. mode couples 
to the bath, such that a fraction of the initial excitation 
energy remains in the system. The same is true of the 
coherences, i.e. the off-diagonal elements in the density 
matrix. 

High initial excitation - These generic features can be 
analyzed in more detail for the case of a high initial ex- 
citation energy (no — n F 3> 1). Then, the expression for 
p can be given explicitly in a relatively simple form: 

o ,(t) — * f«,>»K-»»)'i , , x 

n ,n' 

[pdocay(^o -n,t) + ph ca t(n - n F - 1, n' - np - 1{J$) 
The decay of the excitation is described by 

Pdec ay (m,t) = ^-(u(t) + u*(t)) m e v ^+ v W, (13) 

where m = uq — n = n' — n' may be interpreted as 
the net number of quanta transferred to the bath (Fig. 
13 left). At short times, R < 1, the nonvanishing en- 
tries are Pdecay(M) ~ rt and pdccay(0,i) w 1 — rt, i.e. 



Figure 3: Generic time evolution of density matrix in the limit 
of high initial excitation (Eqs. Q2j~CU> plot for lo' — loq). 



Golden Rule behaviour is recovered (both for relaxation, 
rio = n' Q , and dephasing, uq ^= n' ). In the long-time 
limit we get a stationary distribution, Pdecay^, i) — > 
(2 rn /ml)e~ 2 . 

"Heating" around the Fermi level is encoded in (see 
Fig. El 



p heat (n,n',t) = v(t) n - n ' 



2(mi+ra 2 ) 



(-1) 



m2+m 2 



(14) 

where the triple sum runs over fh\ = max(0, n' + l) . . . oo, 
fh% = . . . oo, mi = max(0, fhi + n + 1) . . . n — n' + fn\ + 
fri2, and we have — fh\ + rfi2 — rri\ + n — n'. In 
the short-time limit, /?hcat {n,n 1 1) approximates to 1 for 
n < -1, to 1 - \v(t)\ 2 for n = -1, to \v(t)\ 2 for n = 0, 
and for n > (to 0(\v\ 2 )), describing the unperturbed 
Fermi sea and the onset of heating. Comparing to the 
full results (Fig. EJ), we find that the limiting case H12|) 
is a very good approximation even for small excitation 
energies. 

Generic coupling - Up to now, we have considered a 
coupling where the particle coordinates enter linearly, 
and consequently only the cm. mode is damped. 

We now extend our analysis to a more generic situa- 
tion, replacing the interaction in Eq. J3J) by: 



N 



2mujt 



(15) 



Now the bath induces transitions between levels n+q and 
n, with an arbitrary (real- valued) amplitude oc f q (which, 
however, must not depend on n). For f± = 1, f g = (q > 
1) we recover the original model. For fq^O all the bo- 
son modes are damped and couple to each other via the 
bath. Formally, the correlators (b q '(t)b^ q J can be written 

in terms of the resolvent of the classical problem of bo- 
son oscillators coupled to bath oscillators (0], compare 

E3)- 

The correlator ($(x', t)(f>(x, 0)^ now contains contri- 
butions for all pairs q,q'- The evaluation of the Green's 
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Figure 4: Decay of populations for a bath inducing transitions 
between levels that are arbitrarily far apart (see text). 



functions proceeds as before. Unfortunately, one has to 
deal with far more terms. However, interesting behaviour 
is already found in the weak-coupling limit, which here 
implies neglecting the effective coupling between boson 
modes that has been induced by the bath, and describ- 
ing the correlator of each boson mode separately as a 
damped oscillation. 

For the case of constant f q = l (up to some cutoff) and 
an Ohmic bath spectrum, the boson correlator decay rate 

(oc qlPP) , see Eq. (JTSJ) equals q 2 T/2. This fits 

the expectation about Pauli blocking: The decay of a 
particle from state np + Sn + 1 is due to transitions by 
1 to 8n levels, and adding up their rates (which grow 
linearly) leads to a total rate oc 8n 2 , consistent with the 
decay rate of the highest boson mode q = Sn that is 
excited by adding this particle. The actual evolution of 
the Green's function is a superposition of decays, with 
rates up to this value. 

An example for the resulting time-evolution is shown in 
Fig. 21 Starting from a state where a single extra particle 
has been added in level no above the Fermi sea, one can 
observe the evolution of the populations p nn (t) (Ohmic 
bath, T = 0). At intermediate times, heating around 
the Fermi level takes place (barely visible, in contrast to 
Fig. |2J). In contrast to the previous case, the relaxation 
towards the N + 1-particle ground state is complete, the 
system is ergodic. 

Conclusions - We have analyzed a many-fermion gen- 
eralization of the single particle in a damped harmonic 
oscillator, illustrating relaxation and dephasing in a dis- 
sipative many-particle system. Using the method of 
bosonization (in the limit of large particle number), we 
have derived exact expressions for the Green's functions 
and discussed them in limiting cases. We have analyzed 
the decay of an excited state created by adding one par- 
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ticle above the Fermi level, where one can observe the 
"heating" around the Fermi level (due to the effective 
interaction between particles), as well as the incomplete 
decay of the excited particle. Finally, we have extended 
our analysis to a more generic type of coupling between 
system and bath, where the system becomes fully er- 
godic. 
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